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The asymptotic behavior as t — > oo and r ^ oo of the 
hard-thermal-loop propagator *D'^'^{t,f) is computed in the 
Coulomb gauge. The asymptotic falloff is always a power law 
though generally different in the deep time- like and space-like 
regions. The contributions of quasiparticle poles and Landau 
branch cuts are computed. The most difficult calculation is 
the contribution of the branch cut in the transverse propaga- 
tor *D'-'(i,r^. For QED this produces a leading behavior of 
order T/r in both the time-like and space-like regions. The 
inclusion of a magnetic mass so as to describe QCD makes the 
leading behavior l/(Tr^), thus improving the infrared conver- 
gence. The asymptotic space-like behavior of all contributions 
(longitudinal and transverse, poles and cuts) is confirmed by 
also computing in the Euclidean formalism and analytically 
continuing. The results are compared will those for free gauge 
bosons at finite temperature. 
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I. INTRODUCTION 

The hard-thermal-loop approximation to high temper- 
ature gauge theories developed by Braaten and Pisarski 
provides a consistent perturbative method for calculating 
processes in which one or more of the external momenta 
are of order gT ^-^ . An excellent review of the entire 
subject and its applications is given by Le Bellac [Q . 

For an SU(N) gauge theory with N^ flavors of massless 
quarks in the fundamental representation, in the high- 
temperature deconfined phase there is an effective ther- 
mal mass for the gluon 



2 fnr Nf\g^T^ 



(1.1) 



For QED at high temperatures vrT 3> vrif^, there is an 
effective thermal mass for the photon 



er/3. 



(1.2) 



However, these effective masses do not screen the electric 
and magnetic fields produced by all time-dependent cur- 
rents. The situation is much more subtle. The phenom- 
ena of dynamical screening applies only to currents with 
frequencies below rUg g. When the frequency satisfies 
fco < rrig then the momentum-space propagator has poles 
at complex values of the wave vector k. The imaginary 
part of k produces the screening. In a different frequency 
regime, this is the reason that AM radio signals are expo- 
nentially attenuated in the earth's ionosphere H. In the 



hard-thcrmal-loop approximation static magnetic fields 
with fco = arc not screened. Only if < /co < "^-g does 
the pole contribution to the hard-thermal-loop propaga- 
tor, *!)'''' P°''' (fco, r), fall exponentially at large r. How- 
ever, there are also branch cuts in the momentum-space 
propagator. How the cut contribution, *Z3^'''^"*(fco, r), 
behaves at large r has not been investigated. 

Recent calculations of Boyanovsky, de Vega, et al have 
provided complementary information about the asymp- 
totic behavior of the hard-thermal-loop propa gator at 
large time and fixed fixed wave- vector k [ 10|Jll| . Those 
authors were interested in non-equilibrium effects and 
found it convenient to use the retarded propagator, 
*D'^ {t, k). Being the Fourier transform of the fixed fre- 
quency propagator, it necessarily contains the effects of 
all frequencies. They discovered that the Landau branch 
cuts are extremely important for the asymptotic behavior 
of *D'^{t, k) and found power-law falloff rather than ex- 
ponential. Similar results for current-current correlators 
were obtained by by Arnold and Yaffe |l2| . 

This paper will examine the asymptotic behavior in 
both time and space of the hard-thermal-loop propaga- 
tor *D^'^{t,r). The propagator considered will be in the 
Coulomb gauge and it will be the time-ordered propaga- 
tor, not the retarded. The asymptotic behavior in the 
limit t ^ oo, r —> oo with r/t fixed is controlled by the 
region in which the energy ko and momentum k are both 
of order rrig, which is precisely the region in which the 
momentum-space propagator differs significantly from 
the free propagator. Reaching the asymptotic regime will 
generally require t ^ ^/f^g and r 3> l/rUg. The region 
far from the space-time origin with r/t < 1 will be called 
the deep time-like region and that with r/t > 1 will be 
called the deep space-like region. 

It is difficult to anticipate the asymptotic behavior. 
The case of a free, massive vector boson is an interesting 
comparison. In the deep time-like region that propagator 
falls like a power, 1/i^/^; in the deep space-like region it 
falls exponentially, (T/r) cxp(— Mr). 

For the hard-thermal-loop propagator the one might 
expect the quasiparticle poles to produce an asymptotic 
dependence like that of a massive boson. In the the deep 
time- like region this is correct, i.e. 1/t^/^. However, in 
the deep space-like region the quasiparticle pole contribu- 
tions to both the longitudinal and transverse propagators 
produces power-law fall-off. The failure of the quasiparti- 
cle poles to produce an exponential fall-off confirms that 
the electric and magnetic fields produced by localized, 
time-dependent currents are not generally screened since 
those sources can contain frequencies higher than rrig. 



It is even more difficult to anticipate the effect of the 
branch cuts in the hard thermal loop propagator since 
free propagators can offer no guide. The most surpris- 
ing result is the effect produced by the branch cut in 
the transverse propagator. The asymptotic behavior in 
either the deep space- like or time- like regions is the same: 



*2?*J™*(i,r) 



-iT 



(^gij + f^f.j^ + 



(1.3) 



In the deep space-like region the T/r behavior is not sur- 
prising since that is the leading behavior of free ther- 
mal propagators. However in the deep time-like region 
the T/r behavior is very non-trivial. In contrast, free 
propagators in the deep time-like region fall cxponcn- 
tially|l||l|. ^ ^ 

In QCD the inclusion of a magnetic screening mass will 
change the above effect of the branch cut in the transverse 
propagator. This change will be c alcu lated. In QED 
there is no magnetic mass and Eq. (L3) applies. 

The paper develops systematically the mathematical 
machinery necessary for the analysis. There are, of ne- 
cessity, many contributions to evaluate. There are two 
propagators, longitudinal and transverse. Each propaga- 
tor has a pole and a cut to consider. Each contribution 
must be evaluated in the time-like and in the space-like 
limits. The emphasis throughout is on extracting the 
asymptotic power-law behavior of the propagators. Ex- 
ponentially small corrections are always omitted. Sec- 
tion II examines the simple example of the free photon 
propagator in Feynman gauge. Section III summarizes 
how gauge boson propagators, free or interacting, are ex- 
pressed in terms of spectral functions. The real subject of 
the paper, the asymptotic behavior of the hard-thermal- 
loop propagator, is confronted in Sees. IV and V which 
treat the longitudinal and the transverse propagators, re- 
spectively. Section VI provides a detailed summary of the 
results. For many readers it may be best to read the re- 
sults in Sec. VI before reading how they were calculated. 

There are four appendices. Appendix A displays the 
hard-thermal-loop spectral functions in other gauges. 
Appendix B proves the lemma that is used in the text 
to determine the asymptotic behavior in the deep space- 
like region. Appendix C computes the asymptotic be- 
havior of *'D^^'^{t,r) in the deep space-like region by us- 
ing the Euclidean time formalism. The sum of the static 
and non-static contributions agree exactly with the sum 
of the pole and cut contributions that are computed di- 
rectly in Minkowski space in Sees. IV and V. Appendix 
D contains details of the most difficult calculation, the 
asymptotic behavior of *!?*•' '^"'(t, r) in the deep time-like 
region. 



II. SIMPLE EXAMPLE 

Before confronting the problem of hard-thermal-loop 
propagators it is useful to examine the much simpler 



problem of thermal propagators for free gauge bosons, 
i.e. those without hard-thermal-loop resummation. The 
reason for this example is to illustrate two points that 
will be important later. However, the free thermal prop- 
agators are important in their own right in two differ- 
ent contexts: (i) In the low temperature regime of QED, 
■kT <^ rUe, no resummation is necessary for the photon 
propagator, (ii) In the high temperature regime of QED 
or QCD the ghost propagators do not require resumma- 
tion and the free thermal form applies H. 

The free propagator is simple to compute in the Feyn- 
man gauge and provides two useful comparisons with the 
hard-thermal-loop results that will come later. In the Eu- 
clidean time formalism the Feynman gauge propagator is 



-D^ir^r-) 



g^^tT 



71^ — OQ 



d?k e**^-*"-*""^ 



{2n) 



(2.1) 



where ojn — iirriT and —(3 < t < (3. The poles at fc = 
dziLUn allow the momentum integration to be done by 
Cauchy's theorem: 



^'Eir,r)^g^^§-^ 



2 y cos{ujnT)e 



The time-ordered propagator in real time is obtained by 
the replacement r — > it: 



2?A'>'(t,f)=g^'' 



iT 
Anr 



1 + 2 yj cosh(ijj„i)e 



n=l 



In the space- like region r>t each term in the sum is expo- 
nentially small. Therefore the asymptotic behavior in the 
deep space-like region is iT/Airr. However, in the time- 
like region t>r the terms in the series grow exponentially 
with n. To evaluate the behavior in the time-like region 
it is necessary to first sum the series in the space-like re- 
gion (the region of convergence) and then continue that 
finite result to the time-like domain. Since the series is 
geometrical, the summation is elementary and gives 



V^^''{t,r)^g^'^-^ 



iT 
47rr 
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1 



1 



p27rT(r-i) _ I ' g2irT(r+t) __ I 



When r — > cx) this still behaves as l/r plus exponentially 
small corrections. The surprise is that when t — > cx) this 
is exponentially small. More generally, for r — > cx) and 
t —f oo with a fixed ratio r/t < 1, the propagator is 
exponentially small. 

There are two lessons to be drawn, (i) The behav- 
ior in the deep space-like region can be easily calculated 
using the Euclidean formalism, (ii) The behavior in the 
deep time-like region can be quite different from the deep 
space-like. Unless one is able to perform the sum over fre- 
quencies, the Euclidean propagator cannot provide infor- 
mation about the asymptotic behavior of the Minkowski 
propagator in the deep time-like region. 



This example is misleading in one respect. For the 
free propagator the leading space-like behavior, T/r^ 
comes entirely from the static n = mode. For the 
hard-thermal-loop propagators the leading behavior of 
*I?°'^(i,r) will come from the non-static modes. 



III. GENERAL STRUCTURE 

Before getting to the HTL propagators, it is useful to 
display the underlying structure of the gauge boson prop- 
agator that holds whatever the approximation. The basic 
two-point functions or thermal Wightman functions are 



^r(^) 



v':^{x) 



-i^T[QA^'{x)A''{{})\ 
-i^T[QA''{f^)A^'{x)\, 



where g — exp(— /3iJ)/Tr[exp(— /3i/)] is the density oper- 
ator at temperature T = 1//3. 

Time- ordered propagator. The propagator that will 
be calculated later is the time-ordered propagator: 

P'^''(i,f) =0(t)I?^"(i,f) + 0(-t)I?^"(t,f). (3.1) 

What is actually known from hard-thermal-loop calcu- 
lations is the propagator in momentum space or equiva- 
lently the spectral function since the time-ordered prop- 
agator can be expressed as 



i^-(fco,fc)=/ ^^^>'^^) 



, 27r ko — (T + ie 



^C^. (3.2) 



o/3fco 



1 



Knowledge of the momentum-space propagator immedi- 
ately yields the spectral function since 



lmD'"'{ko,k) ^-pt"'{ko,k) 



1 



o/3feo _ I 



The Fourier transform of Eq. ( |3.1| ) implies that the two- 
point functions in momentum space are 



D'iriko,k)^~ip'^''iko,k) 



i?f(fco,fc) = -*P^''(fco,fc) 



1 



1 



o/3/co _ 1 



,l3ko 



1 



Knowing the spectral function allows one to calculate the 
two-point functions in space-time by Fourier transform- 
ing: 



v'tr{t,f) = -i 



'D'"<{t,f) 






(2^)4 

(2^ 



(3.3) 



ik-r—ikot 



p'"'{ko,k) 



The asymptotic behavior of this Fourier transform is the 
subject of this paper. For the asymptotic behavior the 
time will always be positive (and large) so that the time- 
ordered propagator is the same as the thermal Wightman 
function. 



Retarded propagator. For certain applications the re- 
tarded propagator is of interest [^ 11 1 . In space-time the 
retarded propagator is 

V^^ (t, r) = 0{t) [Vtr (t, f) - Pii" {t, r)] . (3.4) 

In momentum space this becomes 

da p'"'ia,k) 



^r(^o,fc)= ^ 



, 27r fcn 



a + le 



Knowing the spectral function allows computation of the 
space-time dependence by 



I?^-^(i,f) = -z 



(2^ 



^ik-f—ikot i^LU 



nko,k). (3.5) 



Comparison of Eq. (3.3) with Eq. ( |3.5D shows that the 
time-ordered propagator contains the Bose-Einstein func- 
tion but the retarded propagator does not. As t — > oo 
the time-ordered propagator is more sensitive to the low 
frequency part of the spectral function. 

Imaginary time propagator. The imaginary-time, or 
Euclidean, propagator already used in Sec. II is ordered 
in the variable r: 



n^{r,^ 



(r)I?(:''(-zr,rO +0(-r)I?f (-jr,r). 



where —(3 < t < j3. If one computes the Euclidean 
propagator for < r < /3 and then analytically con- 
tinues T to the imaginary axis (r — > it), the result will 
be 'D'i!^ [t,r). Thus analytic continuations in time relates 
the Euclidean propagator to the time-ordered propagator 
as used in sample calculation of Sec. II. By contrast, if 
one knows the Euclidean propagator as a function of fre- 
quency, D^ {bJn , 'r) , then the continuation cij„ — > ik^ — e 
will give the retarded propagator D^[kQ, r). 

HTL spectral functions in Coulomb gauge. Calcula- 
tions will be done in the Coulomb gauge, which is used 
extensively in applications of hard-thermal-loops H . In 
this gauge the spectral function has the tensor structure 



p°°(fco,fc) = pKfco,A;) 
p'^{ko,k)^{5'^ ~Uy)pt{ko,k), 



(3.6) 



with six components vanishing: p^^ = p^'^ = 0. The lon- 
gitudinal and transverse spectral functions, introduced 
by Pisarski |l[], both have the form 



27r 



Psiko,k) = Zs[6{ko -ujs) - S{ko +ujs)] 
+e{k^~kl)l3s{ko,k), 



(3.7) 



where s denotes either i or t. Both p^ and pt are odd 
functions of kg. The Dirac delta functions come from the 
pole in the momentum-space propagator at the quasi- 
particle energies ztujs{k). The residues of those poles are 
the real functions Zs{k). The continuum contribution fig 
come from the discontinuity across the Landau branch 
cut in the momentum space propagators. 



IV. LONGITUDINAL PROPAGATOR: *r>°°(a;) 

The main subject is the computation of the asymptotic 
behavior of the hard-thermal-loop, time-ordered prop- 
agator *'D^^{t,r) which for positive t is the same as 
the Wightman function *'D'i^ {t^r). This section treats 
*I?™(i, 7^), the Wightman function for gauge bosons with 
time-like polarizations in the Coulomb gauge. In momen- 
tum space it is determined by the longitudinal spectral 
function: 



D^^iko.k) 



Pe{ko,k) 
1 _ e-/3fco ' 



(4.1) 



The spectral function is 



1 



— Pi{kQ,k) = Ze[5{kQ - ujg) - d{kQ + cui)] 
+e{k'~kl)p,{ko,k), 
where the residue function is [y 



Z( = 



k-^iiml-wj + k^y 



and nig is the thermal gluon mass. The longitudinal spec- 
tral function for the cut is m 



(3e{ko,k) = 



3m2 Ni{x, k) ' 



(4.2) 



where x — ko/k and the denominator is 



Niix.,k) 



r 2P 



-2-x In 



1 + x 
l~x 



+ [ttx]\ (4.3) 



The longitudinal spectral function satisfies various sum 
rules PM| . The one that will be useful subsequently is 



' „ l3e{ko,k) 

dko 

-fc ko 



3ml 



k^{3ml + k^) 



2Z, 



(4.4) 



To compute the Fourier transform of Eq. ( [l.l| ) it is effi- 
cient to organize the calculation into two steps: 



Fi{t,k) 



fk_ 
dko e 



^ — iknt 



2tt 1 



-I3h 



■ pe{ko,k). 



The integral over the solid angles of k can be done di- 
rectly: 



-I 
27rV 



*X>>° it,r) = ^^ I dk k sm{kr)Fe {t,k). 



(4.5) 



The function Fg is a sum of the pole and cut contributions 

Fe = F/°''= + F™S where 



Fl°'\t,k)^Z, 
Fr\t,k)~- 



-iuj^t 



giuit 



+ 






(4.6) 



dko 



1 _ e-/3fco 



/3,(fco,A:). (4.7) 



The Wightman function at fixed three-momentum is 
^''^{t,k) = -iFi{t,k). 



A. Asymptotic time-like behavior of *X'2?(a;) 

The deep time-like region requires i — > oo and r — > oo 
with a fixed ratio t/r > 1. 

1. Pole contribution. 
The contribution of the quasiparticle pole is 



+ .T-jOO pole 



it,r) 



2Tr^r 



xZf 



dk k sin kr 



-iujf^t 



^iuj£t 



Extending the integration to negative k gives 



*V' 



,00 pole 



it,r) 



-1 



Att^j 



dk 



xkZt 



^ikr—iutt 



~ikr+iuJtt 



\ — Q-Pi^t \ — g0i^t 



The asymptotic behavior of this will be computed by the 
method of stationary phase [l^jlq] . Define 



(4.8) 



(t>{k) — kr ~ LiJe{k)t. 



At large values of r and t the functions e^*"^ oscillate 
rapidly except at special values of k. Let k be the value 
of the wave vector at which the phase is stationary: 



= 



dcj) 
dk 



duji 



dk 



This requires that the group velocity of the quasiparticle 
satisfy 



duji 
'dk 



t 



< 1. 



(4.9) 



The group velocity of the longitudinal dispersion rela- 
tion is zero at fc = and grows monotonically with k, 
approaching 1 as fc — > oo . Therefore there is always a 
real k which satisfies this. Near the stationary point the 
phase is 



<j>ik) ^ cj^ik) - ^{k - k)H 



d^LOl 



dk^ 



The acceleration (Poje/dk'^ is always positive. To evalu- 
ate the contribution of the stationary phase point, the 
k contour must be deformed into the complex plane so 
that it passes through the real value of A: on a path of 
steepest descent. For the first integrand, ex-p{+i(f>), this 
requires k = k + e^'^'^''^s where s is real and gives a simple 
Gaussian integral: 



dk e* 



g^0-»V4 dse''''""/^. 



For the second integrand, exp(— i 

A: = fc4-pW4 

gives 



the correct path is 



e'"' 's. Performing the two Gaussian integrations 



*I?°"P°'°(t,r) 



-k Zi 



r(27r)3/2 y^ 



(4.10) 



^ — i04-*7r/4 



1 



-Piui 



oP^l 



1 



Since r/t is fixed, the value of k does not change while 
t —r oo. The dependence on the size of t arises only from 
the power behavior l/{rt^^'^) and from the phase factor 
(j) = kr — uji{k)t. 



2. Cut contribution 

To evaluate the cut contribution in the deep time-like 
region it is convenient to replace the integration over fco 
by an integration over the dimensionless ratio x — k^/k 
so that 



*I?S,°™*(i,r) 



Ii{x,k) 



-1 



(2^)2 



dk da;e*'=('-^*)/f(a;,fc) 



k'^(3i{kx, k) 



(4.11) 



The phase 4'(x, k) = k{r — xt) has vanishing first deriva- 
tives with respect to x and k at the stationary point 
X = r/t, k = 0. In order to perform a Gaussian integra- 
tion about this stationary point it is convenient to change 
to new variables u and 9 defined by 



r 
X = - 



1 



t {nigty/^ 



je-in/4 



(?) 



1/2 



ue 



-i0 — iTT/4 



(4.12a) 
(4.12b) 



where u is positive and 6 is between and 2tt. We 
are concerned with the limit rrigt — > oo. As long as 
u ^ {rugty^'^ then the variables x and k are only slightly 
extended into the complex plane. With this substitution 
the complex phase factor becomes a Gaussian: 

^ik(r-xt) _ g-"^ 

This is exact. The Jacobian of the transformation gives 
dkdx = udud9{2/t). The cut contribution becomes 



*V' 



,00 cut 



it,r) 



-1 



27r2rt 



27r 



udu I d9 e '^ h{x,k). 



'0 "'0 

Because of the Gaussian, values of u that are much larger 
than 1 make negligible contribution so that the upper 
limit Mmax can effect ively be taken as cx). The variable 
change in Eq. ( 4.12| ) guarantees that x ^ x — r/t and 
fc — > in the limit rugt 3> 1. (Note that the func- 
tion Ni{x,k) contains momentum via the ratio k/rrig '^ 
{mgty/^ ^ 0.) 

The leading term in the integration would normally 
come from the value of Ii{x,k) at the stationary point. 
However Ii{x, k) vanishes at fc = 0. To find the first non- 
vanishing contribution it is necessary to expand h{x, k) 
in a double Taylor series: 



h{x,k) = 



EE 

7n— n— 



{x- 



'(fc)" 



1 Ir7 ' 



d 



,m-{-n ] 



dx'^-dk-^ 



The same type of series expansion will be used in Sec. V 
A for the transverse propagator. It is therefore worth- 
while to analyze the implications of this expansion with- 
out using the fact that for the longitudinal case I(,{x,k) 
vanishes at fc = 0. 

The first general obs ervat ion is that because of the 
variable change in Eq. (4.12), 



-i'^/i\ m+n , 



,._.,.<„.. (HL_)-(_1_) 



The integral over 9 will vanish unless m = n. Therefore 
the only parts of the Taylor series expansion that will 
survive the 9 integration are 



Ii{x,k) 



E 



— f- 

nlnl V 



t 



92"/^ 



9a;"5fc" 



x,k 



The first few terms in the asymptotic series are 

^2ir 

udu I d9e-''''li{x,k) (4.13) 



TTli{x,0) 



in 

T 



d^h 



dxdk 



x,a 



TT r d-^h 
2i^ld^xd^k 



xfi 



with the next correction of order l/t^. 

For the longitudinal propagator the first term vanishes: 
li{x,0) = as mentioned. The necessary second deriva- 
tive is 



dHe 



dxdk 



2T 

3m2 



d 



1 



dxNe{x,0) 



Therefore the asymptotic behavior of the cut contribu- 
tion in the deep time-like region for rrigt 3> 1 is 



*pU 



\t,r) 



iT 



i-KTri^ rt"^ 



d 



1 



dxNi{x,Q) 



(4.14) 



The quantity in square brackets is rather complicated, 
but it depends only on the ratio x — r/t. Thus the 
asymptotic behavior is entirely given by the explicit pow- 
ers, l/(ri2). This is a more rapid f all-of f than l/{rt^''^), 
the pole contribution given in Eq. (4.10). 



B. Asymptotic space-like behavior of '©''"(a;) 

The limit t —> oo and r —> oo with fixed ratio t/r < 1 is 
the deep space-hke region. As is obvious from the anal- 
ysis of the pole and cut contributions in Sec. A, when 
r/t > 1 there will be no points of stationary phase for 
real values of k and fcg . Stationary phase points at com- 
plex values of k or fco will produce terms that fall ex- 
ponentially. The largest effects will instead be endpoint 
contributions that come from the region of small k in the 
integral 



*Vl°{t,r) 



2n^r 



dkksm{kr)Fe{t,k). (4.15) 



The endpoint contributions are evaluated using the 
lemma proven in Appendix H, which depends upon 
Fi{t, k) being an even function of k that is infinitely dif- 
ferentiable on the the real k axis. Using repeated integra- 
tion by parts it is possible to show that if F^ (i, k) is finite 
at A; = then the integral falls exponentially as r —^ oo. 
If however Fi{t, k) -^ co/fc^ with cq a constant as fc — > 
then the integral is of order 1/r with exponentially small 
corrections. Specifically 

/•OO 

lim / dkksm{kr) Fe{t,k) = h exp small. (4.16) 

This will be used here and again in Sec. V B. 



1. Pole contribution 



The behavior of Ff° °(i, r) shown in Eq. (4.6) at small 
k depends upon the behavior of the quasiparticle residue 
Zt. In the low momentum limit loi —f nig and Z^ — )■ 
TOg/2fc^ SO that 







7?pole 



FP°"(i,0) 



rng_ 
2A;2 



<t). 



where 



z{t) 



1 



^-firug 



opmg 



1 



Applying the lemma gives the asymptotic behavior 



*P°°P°i'=(i,) 



87rr 



■ c(t) + exp smah. (4.18) 



2. Cut contribution 



For the cut, one needs the small momentum behavior 



ofEq. (4.7) 



fc^O: Fr\t,k)^T / dko , 

-k ko 



Pe{ko,k) 



[l + Oikl)]. 



The value of this integral is given by the sum rule in Eq. 
(4.4). At small momentum, Zi/uje -^ l/{2k'^), which 
makes the right hand side finite at fc = 0. Using the 
low momentum behavior of the longitudinal quasiparticle 



energy, uj^ 



3k^/5 + . . ., in the residue Z^ gives 



for the low momentum limit: 



Fr\t,k) 



4 T 
15 m^ 



0{k^). 



Since this is finite, by the lemma displayed in Eq. (4. If) 
*I?°'^'^"*(i, r) does not behave as 1/r in the deep space- 
like region but instead falls exponentially with r. The 
only part of *2?°°(i, r) that does not fall exponentially is 
the pole contribution already given in Eq. ( [4.18 ). In Ap- 
pendix y the same result is obtained using the Euclidean 
time formulation. 



V. TRANSVERSE PROPAGATOR: "V^^i^x) 

The analysis of the asymptotic behavior of the trans- 
verse components of the HTL propagator parallels that 
of Sec. IV. In the Coulomb gauge the transverse compo- 
nents of the momentum-space Wightman function are 



*D^i{k,, k) ^ ^^ f*■^'''l^l, {6^^ - Uy). 



1 



= -/3feo 



(5.1) 



Because the spectr al in tegral for the momentum-space 
propagator in Eq. ( |3.2| ) has poles in k^ at the quasipar- 
ticle energies ±wt(A:) and a branch cut from k^ = —k to 
ko — k, the transverse spectral function has the structure 

— ptiko, k) = Zt [(5(fco - ujt) - S{ko + uJt)] 
+e{k'~kl)(3t{ko,k). 
The residue function is ||l|] 

ijJt{i~ol - k^) 



Zt 



?,mluj1 - (w? - fc2)2 



(4.17) The explicit form of the spectral function for the cut is 







Ptiko, k) 



4 x{l-x^) 
3ml Nt{x,k) ' 



(5.2) 



where x — k^/k and the denominator is 
4fc2 



Nt{x,k)-- 



3TOg 



{l~x^) + 2x^+x{l-x^) Inj^i^V 

LI — XJ- 



TTX{1 — x'^) 



(5.3) 



The spectral function for the cut satisfies various sum 
rules [g|J7|. The one that will be useful here is 



'^ ,, A(fco,fc) 1 2Zt 
ako — 



(5.4) 



k ko k^ ut 

The Fourier transform of *Dy{ko, k) will be written 
d^k 



V'l{t,r) 



(27r)= 



{5'^ ~k'¥)e'''-^Ft{t,k), (5.5) 



Ftit,k)^ -f 



°° dfco e 



-ikot 



27r 1 - e-'3fco 



pt{ko,k). 



Here Ft denotes the sum of pole and cut contributions, 
Ft = ^t' 



pole , pcut 



F™*, where 



-ipolc 



Fr"(i,A:) = Zt 



^ — iujtt 



piujtt 



\ _ g-/3ii't g/3t^t _ 1 



F,= 



\t,k) = J dko- 



'ikot 



'Pko 



Ptiko,k). 



(5.6) 



(5.7) 



As in the previous section, it will be important to know 
the behavior of these two functions at small momentum. 
As fc ^ the quasiparticle energy uJt — 
the residue Zt -^ l/2mg. This implies 



nig and therefore 



0: Fr {t,0) 



pole / 



1 



2m, 



c{t), 



(5.8) 



with c(i) the same function as in Eq. (4.17). For the cut 
contribution, the sum rule in Eq. (5.4) and the fact that 
Zt — > l/2mg implies that 



0: 



F-'(i,fc)->^-^ + 0(fc2). (5.9) 



As fc ^ the right hand side diverges as 1/fc^, which will 
make the subsequent analysis more difficult. 

Angular integrations and tensor structure: It is con- 
venient to first deal with the tensor structure by focusing 
on the angular integrations. The most general rotation- 
ally covariant tensor is 

*V'^ (t, f) = ( - S^^ + 3fV^) G{t, r) nonumber (5.10) 
+ {5'^ -f'f^)H{t,r). (5.11) 

By contracting Eq. ( |5.5| ) with fifj the angular integra- 
tion integration over the solid angle of k can be performed 
with the result 

1 



G{t,r)^-f,f,*V'i{t,r) 



27r2? 



' ,, d /sinfcr\ ^ , , , 

, ^^a;:(^;-)^*(^'^)• 



(5.l2) 



The contraction of Eq. (5.5) with Sij leads to 

l6,,*V^^it,r) 

-i 1"^ 
— — / dk k sin(fcr) Ft (t,k). 
27T^rJo 



H{t,r) 



(5.13) 



The integrands of Eqs. (5.12) and (5.13) both behave at 
small k like k'^Ft(t,k). The previous observations imply 
that as /c ^ 0, k'^Ff°^\t,k) -> and k'^F^'^\t,k) -^ T. 
Consequently the integrations for G and H are conver- 
gent at small k. By using the relation 



Qj,2 



2 5 ,2 
r or 



/sinfcrx 
\ kr ) 



0, 



it is easy to show that the two coefficient functions are 
related by 



d 



^[r^G(t,r))=r'H{t,r). 



(5.14) 



It will oft en be easier to compute the integral for H{t, r) 
in Eq. (5.13) and then use this relation to compute 
G{t,r). Of course it will not determine any part of G 
that behaves like 1/r^. 

The asymptotic behaviors of G and H will be com- 
puted in this s ection , but the results will not be substi- 
tuted into Eq. (5.11) for examination. That examination 
and discussion will be deferred to Sec. VI and the im- 
patient reader interested in the answers rather than the 
computations is encouraged to turn directly to Sec. VI. 



A. Asymptotic time-like behavior of *!'*-' (a;) 

The deep time-like region denotes the limit in which 
t^oo,r^ooata fixed ratio satisfying t/r > 1. 



1. Pole contribution 

The asymptotic behavior of the pole contribution to 
the transverse propagator i s co ntrolle d by a point of sta- 
tionary phase. From Eqs. (5^) and ( 5.13 ) the contribu- 
tion of the quasiparticle pole is 



HP°'%t,r) = 



2n^r 



dk k sin kr 



xZt 



cii^tt 



One can extend the integration to negative k so that 



HP°i°(t,r) = 



-1 

An^r 



dk kZt 



pikr — iuJtt 

1 — e^l^^-t 



-ikr+iuft 



ofii^t — \ 



As before, the method of stationary phase ||lj,|l5[ will be 
used to compute the asymptotic value of this integral. 
Define 



4>{k) — kr ~ LJt{k)t. 



(5.15) 



At large values of r and t the functions e^"^ oscillate 
rapidly except at special values of k. Let k be the value 
of the wave vector at which the phase is stationary. This 



requires that the group velocity of the transverse 
particle satisfy 



quasi- 



dk 



r 



(5.16) 



The group velocity of the transverse dispersion relation is 
zero at fc = and grows monotonically with fc, approach- 
ing 1 as fc ^ oo R. Therefore there is always a real k 
which satisfies the stationary phase condition. Near the 
stationary point, the phase is 



1 



rw2 



(j){k) = (j){k) - -{k ~ ky t 



d^ujt 



dk^ 



+ 



The acceleration d^LUt/dk"^ is always positive. The k con- 
tour must be deformed into the complex plane so that it 
passes through the real value of A: on a path of steepest 
descent. For the first integrand, exp(-|-i(?!)), this requires 
k = k-\-e~'^'^'^s where s is real. For the second integrand, 
exp(— 1(?!)), the correct path is A: = fc -I- e^'^/^s. Perform- 
ing the two Gaussian integrations gives the asymptotic 
behavior: 



7JP°'<=(t,r)- 



-kZt 



(5.17) 



(27r)3/2[l-e 



pZ0— z7r/4 



^ — icji+iir / A: 



^~l3ujt 



pPi^t 



1 



Since the ratio r/t is fixed, this falls like t"'^/^. 

The tensor structure of the transverse propagator in 
Eq. ( |5.11| ) depends on both H and G. Using Eq. ( |5.14| ) 
gives for the asymptotic behavior of G^°^°: 



GP°''=(i,r)- 



i Zt 



(5.18) 



(27r)3/2 



oi0~iTr/4 



-i0+i7r/4 



o/3wt 



1 



To confirm that this is the dominant term in G^°^° it is 
necessary to recognize that the most rapid dependence on 
r comes from the phase, (j> = kr—ujt{k)t. The r derivative 
of GP"'*^ comes from three sources: from the dependence 
of the phase 4> on r, from the 1/r^ power dependence, 
and from the implicit dependence of the stationary point 
k on r: 



dr 



d(j) 9GP°i° 



' y^polc 



dk aGp°'<' 

d^ dk 



Since d(j)/dr — k the first term is larger than the second 
at large r. The first term is also larger than the third 
term because dependence of /c on r is very slow: dk/dr — 



l/{tuj'l), which is of order 1/r. This confirms Eq. ( 5.18 ) 



2. Cut contribution 

To analyze the cut contribution of the transverse prop- 
agator in the deep time-like region it is convenient to 
begin with the double integral representation 



7J™*(t,r) 



r 



dk 

,(2^ 



■ k sin kr 



-ikot 



XI dko — 

-fc 1 



-/3feo 



Pt{ko,k), 



which follows from Eqs. (5/7) and ( 5.13| ). As previously 
noted, the integral over fco approaches T/k^ at small k 
and therefore the integration is convergent for small k. 
Changing variables from /cq to the phase velocity x — 
ko/k and replacing isinkr with exp(ifcr) gives 



iJ™*(t,r) 
It{x,k) 



r(27r)2 
k^(3t{xk,k) 



dk 



1 



-pxk 



dxe*'=('-="*)/t(x,/c) 
1 

(5.19) 



This double integral is quite similar to the integral which 
arose for the longitudinal propagator in Eq. (4.11). The 



stationary phase point occurs ai x — r/t, k = 0. The 
contribution of this stationary phase point will be of or- 
der l/r^, which was the final result for the longitudinal 
propagator in Eq. (4.14). 

It is rather surprising that the stationary phase point 
does not give the largest term in the asymptotic expan- 
sion. The subtlety comes ultimately from the fact that 
there is no magnetic screening. There is a hand- waving 
way to get the correct leading behavior. Reference to Eq. 
( ^.9D suggests making the split 

Fr'{t,k)^^ + sFr'{t,k). 

This gives an elementary integration for the leading term: 



i/'^"'(t,r) = 



27r2rjo 
—iT i 



L + sFr'it,k) 



Anr 27r2 



dk ksm{kr 

dfcfcsin(fcr)(5i^™'(t,fc) 



It will turn out that —iT/Anr actually is the correct lead- 
ing term. However this simple calculation does not re- 
ally prove that because one cannot show that the con- 
tribution of the integral over 5F™* is smaller than T/r. 
The problem is that 5i^™' is too awkward to deduce the 
asymptotic behavior. 

A more convincing analysis is to check qualitatively 
that there is a region of integration whose contribution 
is larger than the point of stationary phase. Thus we 
examine the integrand /f of the double integral Eq. (5.1£) 
in the region where x <C 1 and |fc| <C m^: 



It{x,k) 



small X, small k 



AT 
3m2 



[4fc2/3r 



[tts:]" 



This part of the integrand is largest when k"^ /mi is sim- 
ilar in magnitude to x. The large value of /( can be sup- 
pressed by oscillations in the phase factor exp[ifc(r — a;t)]. 
However, at very small k ^ 1/r these oscillations are neg- 



ligible. Thus the dominant region is k 



Ix < 1/r. 



The approximate value of the integrand here is very large: 
/((x, /c) ~ Tvnir^ The effect of this small region on iy^"' 
can be crudely estimated as 



H''''\t,r) 



dk 



dx [Tm 



y] 



-. (5.20) 
r 



At sufficiently large r this dominates the stationary phase 
contribution, T/{m'^r^). Note also that the T/r contri- 
bution is independent of rrig and thus independent of the 
coupling. 

To perform a precise calculation of the leading and 
sub- leading behavior of ff'^"', the essential step is to sep- 
arate out the dominant contribution to the integrand by 
defining 



j(^)(x,A:) = 



4r 



[4fc2/3?- 



[ttx]' 



l'^^Hx,k)=It{x,k)-l''^\x,k). 



(5.21) 
(5.22) 



Obviously I^^' contains the most important behavior. 
With this split, i7™' = 7J(1) + 7J(2) .^^here for n = 1 
or 2 the integrals are given by 



i?(")(i. 



"(277)2 



dk 



^^g*fc(r-xt)/(n)(a.^fc). 



(1) The dominant contribution will come from H^^' . 
This is computed by direct integration in Appendix O. 
That calculation gives the asymptotic behavior 



if(i)(i. 



-iT 



AiT 



Airr 



(5.23) 



with power-law corrections. Using this in Eq. ( 5.14 ) 
gives the asymptotic behavior of G^^-' as 



G(i)(i,r) 



-iT 
87rr 



AiT 



STT^rri'ir^x 



(5.24) 



(2) The behavior of iJ'^^) j^ ^i^q deep time-like region 
can be obtained from the same type of stationary-phase 
method used for *V^^{t,r) in Sec. IV A. The change 
of variables in Eq. ( 4.12| ) gives 



i7(2)(i, 



rt{2n) 



2tt 



udu / dOe-'' /(2)(x,A:). 



At the stationary point, x — r /t and A: = 0, the integrand 
7^2) ( ^^0) vanishes. The leading term is obtained using 
Eq. Km-. 



H^^\t,r)- 



rt(27r) 



92/(2 



dxdk 



x.O 



Working this out explicitly gives 



i?(2)(i^ 



i2T 



rt'^iiim? dx 



1 



1 



(5.25) 



Nt{x,Q) ^2^.2 _ 
Because of the subtracted term this is finite even aXx — 



0. From Eq. (5.14) the other coefficient G^"^' must have 
the asymptotic behavior 



where fix) satisfies 

df{x) ^ 
dx 



2iT 



SiTm'ir^ 



/(S), 



(5.26) 



dx 



1 



Nt{x,0) 



Integrating this is not automatic because one does not 
have an initial condition. As a; -^ the right hand side 
is of order x"^. Thus f{x) -> /(O) + 0{x^) asx^O. It is 
not obvious that /(O) should vanish. (This would make 
G'2' finite at r = 0, but there is no reason to demand this 
since G^^-* is not finite at r = 0.) Appendix O performs 
a stationary phase calculation of G*^2) with the result 



/(^) 



dyy-r- 

dy 



1 



i-y' 

Ntiy,0) (7ry)2 



In the calculation of Appendix ^ the integration variable 
over y is a remnant of the angular integration. Specifi- 
cally, y = xcos9 where k ■ f = cos 9. Since the quantity 
in square brackets is finite at y = 0, one can integrate by 
parts to obtain 



fi^) 



x{l- 



Ntix,0) n^x Jo 



dy 



Nt{y,0) {iryy 



In this form the entire integral of the last term l/(7rj/)2 
appears important. Because the expression in square 
brackets is finite at ?/ = 0, one can replace the lower 
limit of zero by an infinitesimal value j/o and take the 
limit ya —I- after performing the integration. Doing this 
gives 



/(S) 



x{l 



Ntix,0) 



lim 



vo Nt{y,0) 7r2j/o 



In this form the divergence from the lower limit of the 
integration is canceled by the term —l/{n'^yQ). 

Total cut contribution. In the last form for f{x) the 
term — 2/(7r2a;) when substituted in to Eg . (5.26|) exactly 



cancels the order 1/r^ term in Eq. ( 5.24 ). Consequently 
the sum of G(i) and G(2) is 



G™'(i,r)- 

._^ ^(1 



-iT 



2iT 



Sirr 

-^2) 



STTUlil 



■9{x) 



(5.27) 



Nt{x,Q) 



lim 



^.yf^My.O) 



71'^J/O 



The surviving terms of order l/r^ come entirely from the 
stationary phase point. Similarly when TJ'^' and H^'^' are 
added together the order l/r^ te rm fr om Eq. ( ^.23 ) is 
canceled by the last term in Eq. (5.25) to give 



7J™'(^,0- 



i2T d 



1 



^2 1 



Ntix,0) 



The asymptotic behavior of G"^"' is 



(5.30) 



iJ™'(^,0 



-iT 



i2T d 



Anr rt'^STTm-i dx 



l-x' 



Ntix.Q) 



(5.28) G'^^{t,r)^ 



The corrections to 1/r come from the stationary phase 
point. 

Inclusion of a magnetic mass. It is widely believed 
that non-abelian gauge theories at high temperature gen- 
erate a magnetic mass of order g^T defined by the limit 

u^ ^ limnt(0,/c). 

It is not known how a magnetic mass would change the 
transverse spectral function at nonzero values of fco and k. 
The simplest possibility is that the occurrence of k"^ in the 
transverse spectral function is replaced by the sum k- + 
/i^. Thus the denominator of the new spectral function, 
denoted with a tilde, is 



Nt{x,k)-- 



4(fc2+^2) 



Zm?g 



+ 



{l~x^) + 2x^+x{l~x- 

2 



TTX{1 — X ) 



(5.29) 



The function iy^"' becomes 
-1 



H''''\t,r) 



r(27r)2 



dk 



da;e''=('-^*)/i(a;,fc). 



where the magnetic mass now occurs in the denominator 
of the spectral function: 



it{x,k) = 



k^x{l~x^) 



3m2 [l-e-^'''=]Ntix,k) 



Because of the magnetic mass It{x,k) does not diverge 
at X ~ ^Jk/nig -^ 0. 

The asymptotic behavior of if <="* in the deep time-like 
region can be obtained from the stationary-phase m ethod 
used in Sec. Ill A. The change of variables in Eq. (4.12) 
gives 



7i"™*(t,r) 



rt(27r)2 



udu / dO e ^ It{x,k). 
Jo 



2iT 



STTniir^ 



xil-x"^) 



dy^. 



i-y^ 



Nt{x,0) Jo Ntiy,0) 



(5.31) 



Thus in the presence of a magnetic mass the transverse 
propagator falls like 1/t^ in the deep time-like region. 



B. Asymptotic space-like behavior of *'D^-'{x) 

The limit t — > cxd, r — > oo with a fixed ratio t/r < 1 
is the deep space- like region. As observed in Sec. IV B, 
there are no points of stationary phase for real values of k 
and kg. Stationary phase points at complex values of k or 
fco will produce terms that fall exponentially. The largest 
effects are instead endpoint contributions th at com e from 
the region of small fc in the integral of Eq. (5.12): 



G(t,r) 



2'K^r 



*|:(^)^*(^'^)- (^-^^^ 



1. Pole contribution 



Eq. (5^) shows that at zero momentum, F[°°{t,0) = 
c{t)/2mg. The r derivative may be taken outside the 
integral defining G^°^° so that 



GP°^°(t,r) 



d 



27r'^r dr 



' „ sinfcr ^ , , , 
dk^—Ftit,k) 

'^ 



The integral over fc now satisfies the conditions of the 
lemma displayed in Eq. (4.16) and proven in Appendix 



Therefore in the deep space-like limit 



GP°"=(t,r) 



d 



2T:'^r dr 



1 TT C{t) 



r 2 2mn 



-~ic{t) 
ATTnigr^ 



(5.33) 



The corrections to this are exponentially small. Since 
this is a pure 1/r^ behavior, by Eq. (5.14) the asymptotic 
behavior of the other coefficient function is 



iJP°''=(i, 



exp small. 



(5.34) 



At the stationary-phase point, It{x,Q) vanishes and so 
according to Eq. (4.13|) the leading term is 



il™*(i,r) 



rt{2 



2 (~i'^\ 



d^It 



dxdk 



2. Cut contribution 



For the contribution of the cut, the result in Eq. (5.9) 
motivates separating out the dominant behavior by defin- 
ing 



Working this out explicitly gives 



FrHt,r)^j-^ + SFr'{t,k), 
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where SF^'^^^{t,k) is finite at k 



t 

nves 



0. Using this in Eq. 



G"^"'(i,r) 



' d /sinfcr\ T 
2'K^r n dr 



d 



2Tr^r dr 



dk 



kr 
sin/cr 
kr 



A:2 



Both these integrals are convergent at fc = and at 
k = oo. The first integral may be done exactly. For the 
seco nd, th e asymptotic behavior is given by the lemma, 
Eq. (|4.16|) . The resuh is 



G"^"'(i,0 



-iT 



iT 



Sirr 



Anm'ir^ 



exp small. 



(5.35) 



By Eq. (5.14) the other coefficient function has asymp- 
totic behavior 



i?™*(i,r) 



-iT 



47rr 



exp. small. 



(5.36) 



Inclusion of a magnetic mass. The existence of a 
magnetic mass /i^ would change the asymptotic space- 
like behavior of the transverse propagator. The contri- 
bution of the quasiparticle pole does not change appre- 
ciably. However, for the contribution of the branch cut 
the sum-rule in Eq. (|S.4|) changes to 



dk( 



/3t(fco,fc) 

fco 



nt(o,fc) 



2£^ 



as argued by Pisarski g. In the limit of small k this 
gives 







*" ,, Pt{ko,k) 
dfco ■ 



ko 



0{k''). 



Therefore in Eq. ( 5.35| ) the term of order 1/r is absent 
and the coefficient of the 1/r^ changes so that 



G"="'(i,r) 



iT 
Airr^ 



exp. small. (5.37) 



By Eq. (|I|) 



i/™'(t,r) -^ exp small. 



(5.38) 



Among the exponentially small terms is the familiar 
exp(— /ir). However, the dominant behavior is the l/r^ 
term derived here. In Appendix O the same results are 
obtained using the imaginary time approach. 



VI. SUMMARY AND DISCUSSION 

This section reviews the results obtained in the previ- 
ous calculations and compares them with free propaga- 
tors at finite temperature. The summary is self-contained 



TABLE I. Leading asymptotic behavior of gauge boson 
propagators at finite temperature in the Coulomb gauge. 



Type 



Time- like 



Space-like 



Free massless 
HTL *©'"' pole 
HTL *V°° cut 
HTL 'V'^ pole 
HTL "V^^ cut 



exp. small 

l/t^ 

1/r " 



1/r 

1/r 

exp small 

l/r^ 
1/r '^ 



^This changes to 1/i^ if there is a magnetic mass. 
''This changes to 1/r^ if there is a magnetic mass. 



and may be read independently of the previous sections. 
It is organized differently than the text. First comes the 
asymptotic time-like behavior of thermal propagators; 
then the asymptotic space-like behaviors. The propa- 
gators are always in the Coulomb gauge. 

It is perhaps worth repeating that ghost propagators 
do not contain hard-thermal-loops and therefore are not 
resummed |^ . In the HTL approximation the ghost prop- 
agators are just the free thermal propagators. Their 
asymptotic behavior may be found in |13|. 



A. Asymptotic time-like behavior of gauge boson 
propagators 

In the deep time-like region t and r are both large with 
a fixed ratio t/r > 1. For the free propagator t and r must 
be large compared to 1/2ttT. For the HTL propagator, 
they need to be large compared to l/rUg. The leading 
term in the asymptotic time-like behavior is displayed in 
the first column of Table I. 

1. Free, massless gauge boson. In Coulomb gauge the 
scalar potential A^ (x) is instantaneous and consequently 
I?°*'(t, r) ~ 6{t)/ATTr. The exact space-time dependence 
of I?'^(t, r) can be found in [Q. In the deep time-like 
region t > r ^ 1/{2ttT) the results are 



V"°{t,r)^0 

T>^^ {t, r) — > exp. small. 



(6.1a) 
(6.1b) 



Exponential fall-off in the time- like direction is character- 
istic of free massless bosons and was found in the Feyn- 
nian gauge example in Sec. II. The same behavior occurs 
for spin 0, for spin 1 in any gauge, and for spin 2 gravitons 
]l3| . The exponential fall-off can be understood in terms 
of the calculations performed in Sees. IV A and V A. 
The power-law behavior that occurs in the HTL prop- 
agators comes from a point of stationary phase. The 
stationary phase contribution always has exponentially 
small corrections, which have been consistently omitted. 
However, for free, massless bosons since kg — k, there can 
be no value of k for which the phase factor exp(i/c(r — t)) 
is stationary. The entire answer consists of the exponen- 
tially small corrections. 
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2. HTL longitudinal propagator. From Sec. IV A the 
asymptotic time-like behavior of the longitudinal gauge 
propagator for t 3> l/rUg is 



,Polc/ 



poic/, \ 



*'T-\00 cut 



it,r) 



i3 ' 



(6.2a) 
(6.2b) 



where the coefficient functions are 



a^°'°(t,r) = 



-k Zi 



a^Hx) 



-T d 
"iixrHp-Jx. dx 



pj(^ — z7r/4 



1-e 
1 

m{x,o) 



-piiJi 



~i4>+iTT /A 

afSTJc _ ^ 



The function a^"' depends only on the variable x ~ r/t. 
The function a^° ° also depends on the ratio r/t, since the 
stationary point is the value momentum k at which the 
group velocity of the quasiparticle satisfies duji/dk — r/t. 
In addition a^° ° depends on t and r through the value 
of the phase at the stationary point: (p = kr — uoi{k)t. 

3. HTL transverse propagator. From Sec. V A the 
transverse gauge propagator in the deep time-like region 
has the following asymptotic behavior: 



/.pole/, \ 



(6.3a) 



87rr 

.6;™*(s) 



[-5'^ -f 3fV^) 



(6.3b) 



for times t ^ l/"T-g- The hard-thermal- loop propaga- 
tors do not contain a magnetic mass. Therefore Eqs. 
( |6.3|) apply to high temperature QED. For high temper- 
ature QCD the inclusion of a magnetic mass changes the 
asymptotic behavior as discussed in the next paragraph. 
The coefficients for the pole contribution are 



hT\x) 



\x) 



Sirra'pc^ \ yo~*o 



dy- 



1 



y 



1 



y„ Nt{y,0) TT^yo 
x{l-x'^) 
~ Nt(x,0) 



-2T d 



Sirm-ix dx 



1-x^ 



Nt{x,0) 



In Sec. V A the calculation of the leading term, T/r, in 
*X>ij c"t(i^ f) was lengthy. However, it has a profound im- 
plication for the behavior of electrons. As shown by Pis- 
arski the electron propagator dressed with one HTL pho- 
ton has a damping rate that diverges on the electron mass 
shell ^. In space-time this comes from the T/r. Blaizot 
and lancu showed that summing the most infrared- 
sensitive contributions to the electron propagator pro- 
duces a time dependence of the form exp[— aTiln(i)], 
rather than exponential [18[. Boyanovsky, de Vega, et 
al found the same behavior in scalar QED ||l(]| , [ll|] . 

4-. HTL transverse propagator with magnetic mass. 
For QCD the HTL propagators must be modified so as 
to include a magnetic mass /i^. This does not change 
the asymptotic behavior of the pole contribution. As 
shown In Sec. V A it does change the cut contribution 
by eliminating the term of order T/r and modifying the 



coeffi cients of the l/r^ terms. Using Eqs. (5.30) and 



(Plgi 



*P'^™*(i,f) 



' t^ 

.6J'="'(x) 



(-(5*^+3fV^') 

(ju _ f-^fjy (6.4) 



The coefficients depend on the magnetic mass through 
the denomina tor o f the transverse spectral function Nt 
given in Eq. (|^: 



bl^'ix) 



b'r'c^) 



2T 



a 
-2T d 



d-y-F, 



l-y2 x{l-x'^) 



37rm2xH^o '^Nt{y,0) Nt{x,0) 



3nmjx dx[Nt(x,0) 



l~x' 



The more rapid fall-off in Eq. (6.4) than for QED implies 



that the quark damping rates are infrared finite 



,polc 



{t.r). 



iZt 



lfV(27r)3w: 



^i(j)—iiT /A 



'i4>+i7v/A 



alSiUt 



1 



B. Asymptotic space-like behavior of gauge boson 
propagators 



,/polc 



it,r)-. 



-kZt 



?V(2^ 



377?" 



^iip—iir/A 



1 



-pUJt 



. — ifp-ViTT / A 



These separately on t and r through the phase (p. The 
coefficients for the cut contributions depend only on the 
ratio X = r/t: 



The deep space-like region requires t and r both large 
with a fixed ratio of t/r satisfying t/r < 1. For the free 
propagator r and t must be large compared to l/27rr; 
for the HTL propagator they must be large compared 
to l/rUg. The leading term in the asymptotic space- like 
behavior is displayed in the second column of Table I. 
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1. Free, massless gauge boson. In the Coulomb gauge 
I?°*'(t,r) = 5{t)/4:TTr. The asymptotic behavior of the 
transverse propagator in the deep space-hke region can 
be found in |13|] 



I?°"(i,r) 







-iT ^ 



Svrr 
iT 



(S'j +fif^) 



Sirrx 



(6.5a) 



(6.5b) 



There are omitted terms of order 1/r^ and l/r^. 

2. HTL longitudinal propagator. The asymptotic 
space-hke behavior of the longitudinal gauge propagator 
computed in Sec. IV B can be summarized as follows: 



'•X>0"P°i<=(i,r) 



where 



"V' 



c{t) 



,00 cut 



{t,r) 



.mgc{t) 
87rr 
exp. small 



1 



-Pmg 



1 



o/3mg 



(6.6a) 
(6.6b) 

(6.7) 



The cut contribution is exponentially small as r — > cx). In 
Appendix |Q the same results are obtained in the imag- 
inary time formalism. From that viewpoint, the static 
n = mode is exponentially small because of Debye 
screening. The non-static contributions give exactly the 
same asymptotic behavior as the pole contribution in Eq. 

3. HTL transverse propagator. From Sec. IV B the 
asymptotic behavior of the transverse propagator in the 
deep space-like region is given by 



*X'*JP°i'=(i,r) 



-ic{t) 



SvTTTigr-^ 



(-(5*^'-f 3fVJ') 



-iT 

Sirr 
iT 



Anm'^r^ 



[-5'^ +ir'f^). 



(6.8a) 



(6.8b) 



It is rather unexpected that the leading term T /r comes 
entirely from the branch cut in the transverse propa- 
gator. The subleading term l/{Tr^) comes both from 
the branch cut and from the pole. In Appendix both 
the leading and the subleading terms are obtained in the 
imaginary time formalism. In that approach the leading 
T /r term comes from the static n — Q mode. The sub- 
leading contribution comes from summing over the non- 
static modes. This applies to high temperature QED, 
where there is no magnetic mass. 

4-. HTL transverse propagator with a magnetic mass. 
For QCD a magnetic mass /x^ must be included. The 
asymptotic behavior of the pole contribution does not 
change. As shown in Sec. V B, the term of order T/r 
in the cut contribution is absent and the asymptotic be- 
havior becomes 



*P*J'P°i°(t,f) 
*P*^'™'(t,f) 



-ic{t) 
Snmgr^ 



{~d'^ + 3f'P) 



(6.9a) 



iT 
Airr^ 



1 



1 



m2 /^2 



(-(5*^+3fVJ). (6.9b) 



Exactly the same asymptotic behavior is obtained using 
the imaginary time formalism in Appendix O. In that ap- 
proach the l//i^ term comes from the static n = mode 
and the other terms come from the nonstatic modes. 



C. Light-cone behavior 

The asymptotic behaviors computed for t/r > 1 and 
for t/r < 1 do not allow any conclusion about the behav- 
ior of the HTL propagators on the light cone, t/r = 1. 
To deduce the light-cone behavior, consider the inverse 
Fourier integral over space-time: 



*^M^ 



(fco,fc) 



d^xe'^''*V'^''it,r^. 



(6.10) 



After performing the angular integrations, the phase can 
be expressed as 

kot -kr= -{ka + k)(t-r) + -{ko-k){t+r). 

In the limit feg -|- /c ^ oo with kg ~ k fixed the integration 
is dominated by the light cone i — r — > 0. At large energy 
and momentum the HTL propagator reduces to the free 
thermal propagator in the Coulomb gauge 0| . Thus 



lim 

ko-\-k — *oo 



*D'iiko,k) 



.27r£(fco)(5(fcg-fc2) 

1 _ g-pka 



(S'- 



k'k^] 



Since pko — > oo, the limit is actually independent of tem- 
perature. Because the momentum-space limit k^ -{- k —> 
cx) with ko — k fixed is controlled by the light cone 
t— r —f this implies that the HTL propagator in space- 
time should have the same light-cone behavior as the free 
propagator in space-time. Near the light cone the singu- 
lar behavior of the free Coulomb gauge propagator is 



t^r: V'l(t, r) 



■{5'^ -f'f^) 



87r2i 



B7r2r(i 
■ln(t-r)(-^'^+3fV^) + 



with terms that are finite at t — r omitted. Therefore the 
HTL propagator *T>^^ {t, r) in the Coulomb gauge should 
have these same singularities at t = r. 
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APPENDIX A: SPECTRAL FUNCTION IN 
OTHER GAUGES 

All calculations in the paper are done in the Coulomb 
gauge and the spectral function has the structure given in 
Eq. (3.6). In any other gauge the various components of 
p'^" (K) can also be expressed in terms of the longitudinal 
and transverse spectral functions. 

A particularly simple choice is the temporal gauge, 
A'^ = 0. There are seven vanishing components of the 



r,00 



spectral function: p"" = p 
nishing components are 



P3 



nJO 



0. The only nonva- 



l-il-j 

p^3{K) = iS'^ ~ k'y)pt{K) + -^^Pi{K). 

In the Landau gauge the various components of the 
spectral function are 

p'^W = J0MK) 

In more general covariant gauges the spectral function 
has an additional term K'^K'^d'{K^) that must be added 
to the Landau gauge result. 



APPENDIX B: LEMMA ON ASYMPTOTIC 
SPACE-LIKE BEHAVIORS 

In the deep space-like region, r — > cx) with a fixed ra- 
tio r/t > 1, there are no points of stationary phase on 
the real k axis. For the pole contributions the stationary 
phase condition dco/dk — r/t cannot be satisfied for real 
k. For the cut contributions the stationary phase condi- 
tion k^/k = r/t cannot be satisfied since k^ < k oti the 
Landau cut. In the complex k plane there can be points 
of stationary phase but they will give contributions that 
fall exponentially with r. However, the leading asymp- 
totic behavior will typically be be a power-law falloff that 
is not related to the stationary phase points. The generic 
integral to be investigated is 



V{t,r) 



,, sinfcr „ , , , 

dk— J{t,k). 

kr 



(Bl) 



The function /(i, k) will always be an even function of 
k and infinitely differentiable on the real k axis. The 
lemma will show that asymptotic behavior is controlled 
by the endpoint k = Q. Specifically in the deep space-like 
region 



V{t,r) 



2r 



f{t, 0) -I- exp. small. 



(B2) 



This is the form in which the lemma is used in Sec. V B. 
However in Sec. IV B the naturally occurring integral is 
of the form 



V{t,r) 



1 



dkksm{kr)F{t,k). 



By renaming F{t,k) — f{t,k)/k'^ it is clear that the 
asymptotic behavior of the integral will be co7r/2r where 
F{t, k) — > co/fc^ as fc ^- 0. It is this second form that 
occurs naturally in Sec. IV B. 

To isolat e th e l/r contribution, write the original in- 
tegral Eq. (jmj) asV^Vf +Vg where 



Vjit,r) 
Vgit,r) 



dk^^e 
kr 

sin(fcr) 



k 

oc 

dk 







r 



"'/(i,0) 
9{t,k), 



with g{t, k) defined by 



1 



7(i,fc) = -[/(t,fc)-e-"V(t,0)]. 



(B3a) 
(B3b) 

(B4) 



The parameter a exactly cancels in the sum I?/ + TDg. 
(1) The integral Vf is elementary and gives 

l?/(i,r) = /(i,0)itan-i(r/a). 

Since a is a fixed parameter and r — > cx), one can expand 
the inverse tangent to get 



l?/(i,r) = /(i,0)-[f-^ 



{~lf (a 



1=0 



2£ + 1 V r 



2£+l 



(B5) 



(2) To investigate Vg it is useful to use repeated inte- 
gration by parts. First write "Dg as 



Vg{t,r)^-^ 



„ dcos(kr) , ,, 
dk ^1^ — -g{t,k). 



An integration by parts gives 



Vg{t,r) 



1 



.9(^,0) 



1 



dk 



dk cos{kr 



dgjt, k) 
dk ' 



There is no surface term from the upper limit since 
g{t, co) = 0. Next write cos(fcr) = r~^dsm(kr)/dk and 
integrate by parts. Because of sin(fcr) there will be no 
contribution at fc = 0: 
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Vg(t,r) 



lim 



dk sin(fcr 



(Pg{t, k) 
dk^ ' 



Continuing this process gives 



'Dg{t,r)^Y. 



f=0 



Ml! 
(-ir 



2e„-\ 



d^'g 



dk 



21 



k=0 



,v.2n+3 

' -'0 



dk sin(fcr) 



^2«+2 



dk^n+2 ■ 

From the definition of g{t, k) the even derivatives are 



d''g 



dk^'^ 



1 



fc=0 



2£+l 



d^'+^f 
dk^f'+'^ 



21+1 



fim 



fe=0 



The first term vanishes since f{t, k) is an even function 
of k. Therefore 



Vg{t,r) 



«'-»)7t^(7 



£=0 
(_l)n roo 



21+1 



„2n+Z 



^2n+2 



(3) Adding together Vf and Vg results in a cancehation 



of the powers r , r 



-2n-2 



and gives 



2?(t,r)=/(t,0)-+^ ' 



2r r2"+3 7q 

-/M)^ £ ^(7) 



(ifc sin(fcr' 



^2„+2 



9 



dk^n+2 



21+1 



(B7) 



t=n+l 



This formula is exact. The right hand side is independent 
of the choice of a and of n. (For example, one could set 
a = and eliminate the summation.) The first correc- 



tion to 1/r is of order 1/' 



,2n+3 



Since /(i, k) is infinitely 



differentiable, n may be chosen arbitrarily large. This 
proves that the corrections to l/r fall faster than any 
power of l/r, i.e. fall exponentially with r. The asymp- 
totic behavior is therefore as given in Eq. (B2). It is, 
of course, possible for /(i, 0) to vanish. In that case the 
leading behavior falls exponentially with r. 



APPENDIX C: ASYMPTOTIC SPACE-LIKE 
BEHAVIOR USING IMAGINARY TIME 

The analysis of the hard-thermal-loop propagator in 
the deep space-like region showed that both the longi- 
tudinal and transverse propagators fall like l/r. In the 
longitudinal case the coefficient of l/r depends on the 
thermal gluon mass rUg and on time. This contribution 
comes entirely from the longitudinal quasiparticle pole. 
For the transverse propagator the l/r contribution is the 
same as for a free, thermal field theory. Rather surpris- 
ingly, this simple contribution comes entirely from the 
cut in the transverse propagator. 



This Appendix computes the asymptotic space-like be- 
havior of the longitudinal and transverse propagators 
using the imaginary time formalism. The results ob- 
tained are identical with Eqs. ( |6.6[ ), (6.8), and (|6.S|). 
The asymptotic behavior of the longitudinal propagator, 
which comes from the quasiparticle pole in the real-time 
method, will come from all the non-static modes in the 
imaginary-time method. The asymptotic behavior of the 
transverse propagator, which comes from the branch cut 
in the real-time method, will come from entirely from the 
static mode of the imaginary-time method. 



1. Asymptotic space-like beiiavior of *I''"'(i,r) 

In imaginary time the propagator for the A^ compo- 
nent of the vector potential in the Coulomb gauge is 



*I?°?(r,r)=^^^ 



d^k 



gik-f—iLUn'T 



n— — 00 



(2^)3 fc2+^,(^c^„,fc)^ 



(CI) 



where a;„ = 2TmT and —f3<T<(3. The denominator 
depends on ujn only through the self-energy. Separate 
this into the static, n=0, mode and the nonstatic modes 



V%'{T,r) ^VZoir,r) +V:%{r,r) 



(C2) 



The static value, n£(0, k) — 3™^, produces the standard 
Debye screening result: 



Vll,ir,r)^^T 



d^k 



^ik-r 



iT 



r-N/a 



(27r)3 fc2 + 3to2 47rr 



This exponentially small contribution is not the domi- 
nant one. 

The nonstatic terms will turn out to dominate at large 
separations. The reason is that for n 7^ the propagator 
denominator vanishes at small k: 

o 

fc^O: k^ + Ue{iujn,k)^k^(l + —^)+0{k^). 



Therefore the leading nonstatic contribution in the deep 
space-like region is 



Vl%,iT,r)^tT 



d'k e-^-^ — 



{2n) 



E 



n— — 00 '^ y 



The integration gives l/(47rr) and the sum is elementary. 
The sum converges for r real and gives 



K%ir,r) 



Sirr 



1 



-Prrig 



p/Srng 



1 



This can be continued from Euclidean time t to real time 
t by the replacement r —^ it. Thus use the function 



cit) = 



-pmg 



gl3m,g 



(C3) 
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that was introduced in Eq. (4.17). The space- like asymp- 
totic behavior of the Minkowski two-point function is 



l?00(i,r) 



At) 



Snr 



exp. small. 



(C4) 



This coincides with Eq. (6.6). In the real-time calcu- 
lation this result comes entirely from the quasiparticle 
pole. 



2. Asymptotic space-like behavior of *'D^-'{x) 
The transverse propagator in Euclidean time is 



*V%ir,r)^~^TJ2 



n— — oo 



(2^)3 ujl + k^ + ntiiu;n,k)' 



As before, separate out the static, n — mode from the 
others: 

Because the transverse self-energy vanishes at n = the 
static propagator is the same as the free propagator: 



-iT 



Ak-r 



Sirr 



(27r)3 fc2 



■{6'^ -k'k^) 



(C5) 



This will turn out to be the dominant term. It coincides 



with the cut contribution in Eq. (6.8). 

The remaining sum over non-zero modes can be writ- 
ten 

X'^o(^'^1 = (~S''^' + V^V^)W„^o(r,r), 
where the new function is 



d^k 



^ik-r—iuJnT 



W^Mr.r)- ^TJ^J (27r)3 p [c.^ + p + n,(^c.„, fc)] 

The dominant contribution at large r comes from fc ^ 0. 
In this limit nt(iti;„,0) = m^ and thus 

/d^h „ik-r—iuj„T 
TTT-^ t.2 r 2 I — iT 

_ -iT ^ e-*'^"^ 
47rr ^ Loi -I- m}„ 

The sum is convergent for real t and gives 



I^„5^o(T,r)-> 



'iT 



Airr 



-1 1 



771^ 2m„T 



(r 



-fSnig 1 — e^" 



The continuation from Euclidean r to real t requires the 
replacement r -^ it. so that 



WnMt^r) 



iT 
Airr 



c{t) 



2mgT 



Performing the spatial derivatives gives for the asymp- 
totic behavior of the non-static propagator 



*q~,tj 



KLoit,r) 



iT 
47rr3 



c{t) 
2mgT 



i-S'^ +3f'f^). (C6) 



The leading contribution is the static term, Eq. (|C5|). 
This sub-leading term agrees with Eq. (6.8). In the real- 



time calculation the sub-leading term comes both from 
the quasiparticle pole and from the branch cut. 

Inclusion of a magnetic mass. A magnetic mass elim- 
inates the 1/r contribution of the static mode but intro- 
duces a new l/r3 static contribution. To see this, define 



W„=o(r) = ~iT 



d^k 



ik-r 



(27r)3fc2[p + nt(0,fc)] 
The angular integrations give 



Wn=o{r) 



-iT 
'2^ 



dk 



sin kr 



1 



kr k^+nt{0,k)' 



As r — > (X) this integral is dominated by the region k 
One expects that llt{0, k) ^ fi^ + k"^, in the limit k 
The asymptotic behavior is 



0. 
0. 



Wn^oir) 



-iT 



-iT 



dk 



sin kr 1 



kr k^ + fj? 
, exp. small. 



The exponential fall-off is of the form exp(— /ir). Differ- 
entiating this gives for the static part of the transverse 
propagator 



*^lLoir) 



'iT 



i-S'^ + 3f'f^). 



47rr3/i2 
Adding the non-static contribution from Eq. (|Cq) gives 



'V'\t,r) 



iT 
47rr3 



1 

7 



1 



c{t) 
li 2mnT 



(C7) 



This agrees fully with the real-time result in Eq. (|6.9D. 



APPENDIX D: ASYMPTOTIC TIME-LIKE 
BEHAVIOR OF *D!f'^"'(t,f) 



The tensor structure of the transverse Wightman func- 
tion is 

*V% (i, r) = ( - S'^ + 'Sf'P) G{t, r) nonumber (Dl) 
+ {S'^ -f'f^)H{t,r). (D2) 
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The contribution of the Landau branch cut is particularly 
delicate to calculate. In Sec. V A the cut contribution 
was split into a sum of two parts: 

G™*(t,r)-G(i)(t,r)+G(2)(i,r), 
i?™'(t, r) = i?(i) (i, r) + ij(2) (^t, r), 

in which G'"'^^ and iJ^-'^^ were of order l/r plus order 1/r^, 
whereas G'-^-' and H'^^^ were entirely of order 1/r^. This 
appendix contains the actual calculations of H^^' and 
G*^^-* that were not presented in Sec. V A. 



1. Asymptotic behavior of J/'^'(f,r) 

The leading behavior of the cut contribution comes 
from the integral H'^^\ defined in Sec. V A as 



H^^\t,r)^ 



-1 



^(27r)2 



dk / dxe*'=('^~^*)/(i)(a:,fc), 



with I^^^ given in Eq. (5.21). To compute i?'"^^ it is 
convenient to split the integration over x into three parts, 



with A, B, C defined as 

dk / da; e*'=('^'^*) I^^Hx,k), 

-oo Jo 
/•oo /•! 



(D3) 



dk / dxe*'=(''+='*) /(i)(x,/c). 



G 



At x' = r/i the phase changes sign. The k integrations 
will be performed using Cauchy's theorem, since the only 
singularities in k of the integrand are simple poles. The 
location of the four poles are 

where /ci, A:2, fcs, k^ are located in the first, second, third, 
and fourth quadrants, respectively. The function I^^' can 
be written in terms of partial fractions as 
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"[ [^ '^j)'^^ 



2 • 
3 



(A) In integral A, the quantity r — xi is positive so 
that the integration over k can be computed by adding a 
semicircle of infinite radius in the upper-half of the com- 
plex k plane. Evaluating the k integration by Cauchy's 
theorem gives 



. T [^ dx 
A = - — 



2 ./o X 



ik-iir—xt) ik2{r—xt) 



Using the values of ki and fc2 converts this to 
rdx 



A^iT I —e 
/o X 



-mV^ir-xt) gj^ [^^(^ _ 2.i)j ^ 



where m = m-g\/37r/8. For r —> oo and t ^ oo with 
r/t < 1, the integrand falls exponentially except at the 
two endpoints a; = and a; = a;. The contribution of 
these endpoints will be called Aq and A- and will fall 
like a power of r. Thus with exponentially small contri- 
butions omitted, one can say A — Aq + A-^. 
To compute Aq put x = e^z^ where 



1 _ 1 / 8 
ar m„r VStt 






8 \i/2 



ar TTigr VStt/ 
The the contribution of the lower limit is 



'■ dz 



(D4) 



Aq ^2iT — €-'-+' '' /^ sin [z ~ e'^z'^/x] . 



Expanding the integrand to order e^ and integrating gives 



Ao=iT 



2 X 



(D5) 



To compute the contribution of the upper endpoint, 
a; = a;, return to the exact expression for A and change 
variables to a: = x{\ — ez). Then 



A-7 = iTe 



dz 



^-~z^x(\~e7.) 



sm 



z\/ x(\ — ez) 



JO l-ez 
Expanding the integrand to order e^ and integrating gives 



A^r^iT 






(D6) 



(B) In integral B, the quantity r—xt is always negative. 
Therefore the k contour can be closed in the lower half- 
plane and evaluated by Cauchy's theorem. The poles at 
fca and k^ contribute and give 



z /^ X 



iksir-xt) _ pikii^ — xt) 



Substituting the value of the roots k^ and ^4 gives 

B^iT ^f.rr,V^{r-xt) ^j^ U^(r _ xt)] 

Note that r — xt < in this integral. The integrand 
is exponentially suppressed except at the lower endpoint 
X = X. This contribution B- will fall like a power of r. 
To compute this, put x = a;(l -I- ez) so that 



B-r 



'iTe 



dz 



1 + ez 



e-^(i+")''%in[z(l + ez)i/2] 
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Expanding this for small e and integrating gives 

-2 



B-^iT 






(D7) 



(C) In integral C, since r + xt is always positive, the k 
contour can be closed in the upper-half plane and evalu- 
ated with Cauchy' theorem: 



C 



T r'dx 



iki(r+xt) _ ik2{r+xt) 



Using the values of ki and ^2 gives 



C = iT 



g-mv^(r+xt) gjj^ [^^(^ + xt)] , 



where m = nig -^/Stt/S. The integrand falls exponentially 
except at the lower endpoint x — 0. To compute this 
contribution, put x — e^z^ to obtain 



Co = 2iT 



' dz 



g-.-s^.V- sin [^ + eV/x] 



Expanding in e and integrating gives 

-2 



Co-zT 



2 a; 



(D8) 



Total for H^^' . The sum of the endpoint contribu- 
tions from X w is 



Ao -H Co = iT 



n + Ole^ 



The contributions from x p^ x sum to 



A^ + B^^ iT 



2— + 0(e-^) 
. X 



(D9) 



(DIO) 



It is quite likely that the 0{e^) terms actually cancel in 
this sum, but it will not matter below. Using these in 
Eq. (D3) gives for the asymptotic behavior of H^^^: 



H^-^\t,r) 



-iT 



AiT 



Anr Sn^mir^x 



oi^ 



(Dll) 



This value is used in Eq. ( 5.23 ) 



2. Asymptotic behavior of G^^' (i, r) 

In Sec. V A the asymptotic value of ff^^-' was com- 
puted by the method of stationary phase. This did not 
completely determine G*^^^ and so it too must be com- 
puted by the method of stationary phase. 

It is convenient to start with the original definition of 
G in Eq. ( 5.12 ). If the angular integration over z = k ■ f 
is not performed, the cut contribution is 



G™*(i,r) 



167r2 



dzil-z')g'^\t,r,z), (D12) 



where 



g''''\t,r,z)^ dk da;e*'=(''^-^*)fc/t(x,/c). 



with /( defined in Eq. (5.19). If the z integration is 
performed at this stage it will change kit to It/k and 
therefore prevent converge nce at fc = 0. Separating /( — 
/(I) -|- /(2) as given in Eq. ( 5.22 ) gives 



g^'^\t,r,z)^ dk dxe*'=('^^"^*)fc/(2)(a;^/fc). 

J-oo J ~1 

The integrand fc/2 vanishes like k^ at small k. The sta- 
tionary phase point \s x = rz/t, k — 0. Making the 
variable change analogous to Eq. ( [4.12| ) gives 



g^'^\t,r,z)^ - udu de'e~"'/c/(2)(a;,A:). 

t Jo J a 

The Gaussian integration can be performed using Eq. 
(4.13) to obtain the asymptotic behavior 



t;(2)(i,r,z)^^ 



d^[ki(^\x,k)\ 

dx'^dk'^ 



x,k 



This simplifies to 



g''>{t,r,z) 



-8ttT 92 



3m2t3 ^2 



1-x 



Nt{x,0) (^x)2 



(D13) 



The stationary point depends on the angular variable z 
since x = xz w ith x = r/t < 1. When t/(2) jg substituted 
into Eq. (D12) the integration over z must be performed. 
It is convenient to change integration variables from z to 
y — xz so that 



G(2)(i,r) 



167r2 



dy 



(l-|J)e'(i,r,z). 



Using the asymptotic form for t/(2) gives 
iT 



G^^\t,r) 



iiTirrr^ Jo 



dy{x -y) 



dy^ 



l-y2 



Ntiy,0) {Tryy 



(D14) 



The quantity in square brackets is an even function of y 
that is finite at y = 0. The first y derivative vanishes at 
y = 0. Therefore one can integrate by parts to obtain 



G(2)(f, 



2iT 



dyy 



d 



l-y2 



iirm^gr^ Jo dy 

This is the result used in Eq. ( 5.26| ) 



Nt{y,Q) [iryf 
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